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Abstract. Using the density matrix renormalization group technique, we evaluate the low-energy spectrum 
(ground state and first excited states) of the anisotropic antiferromagnetic spin-one-half chain under mag- 
netic fields. We study both homogeneous longitudinal and transversal fields as well as the influence of a 
transversal staggered field on opening of a spin-gap. We find that only a staggered transversal field opens 
a substantial gap. 



PACS. 75.10.Jm Quantized spin models 

75.50.Ee Studies of specific magnetic materials ( Antiferromagnetic ) 



1 Introduction 



1.1 Model and method 



Recently, the properties of low-dimensional quantum spin 
systems in longitudinal, transversal and/or staggered mag- 
netic fields have become of considerable interest. To give 
an example, the two-dimensional compound SrCu2(B03)2 
is a realization of the Shastry-Sutherland model H, close 
to quantum-criticality, with a spin gap of 31 K which 
shows magnetization plateaus in an external field 

Most spin-1/2 systems show little anisotropics in the 
magnetic exchange. The discovery that low-dimensio- 
nal magnetic excitations in the rare-earth compound Yb4 AS3 
can explain the large linear specific heat coefficient 7 in 
this low-carrier half-metal opens the possibility to 

study in deeper detail the properties of a rare-earth quan- 
tum-spin-chain-system with its enhanced magnetic aniso- 
tropies §p. 

Inelastic neutron scattering experiments § on Yb4As3 
found a gap to (all) magnetic excitations opening in the 
presence of an external magnetic field, confirming a pre- 
diction H by Schmidt et al based on an interpretation 
of previous specific heat data ||. Several proposals have 
been made in order to explain this very unusual behavior 
ppf . The first model || is based on inter-chain interac- 
tions. The second model [[ll] is based on the observation 
that a staggered Dzyaloshinsky-Moriya (DM) interaction, 
which generates an effective staggered g-tensor, is allowed 
HQ in the 4/-compound Yb4As3. It is know that a stag- 
gered g-tensor leads to a gap in an external field |l2],|l3|]. 
The third model |l(| , based on a mean-field analyze of the 
anisotropic spin-chain, proposes that a gap opens in the 
presence of a uniform transversal magnetic field. 

Here we will analysis the two latter proposals by a sys- 
tematic DMRG-studies of the relevant models. We find 
that only the effective staggered g-factor-model is able 
to explain the field-dependent opening of a spin-gap in 
Yb 4 As 3 . 



The magnetic properties of Yb4As3, in the absence of an 
external magnetic field, are well described by an antiferro- 
magnetic Heisenberg spin-1/2 chain. Switching on the ex- 
ternal magnetic field, experimental data shows the open- 
ing of a gap in the low energy excitation spectrum. How- 
ever, the standard Heisenberg model in an applied field 
remains gapless from zero magnetic field up to the satu- 
ration magnetization. The anisotropic Heisenberg model 



L-l 



H = ^ J xy (SfSf +1 + S i S i+1 ) + J Z S*S* +1 

L L L 

- E h * s ? + E K ta9 (-iy sf - E fi'S* a) 



i=l 



i=l 



with an uniform transversal field h x , a staggered transver- 
sal field hg tag and a longitudinal field h z incorporates all 
features proposed to be relevant for Yb4As3. The 

magnetic fields appearing in ([]]) include the gyromagnetic 
g-factors. 

The staggered transversal field in (Q) is induced by a 
staggered Dzyaloshinsky-Moriya interaction given by the 
termY^t-l^D-tSi x S»+i). Setting!) = |D| = J z sin(20) 
the DM-term can be eliminated [^2[Q by a rotation around 
D by an angle 9 leading to h x tag = s'm(9)h z , which can 
be interpreted as an effective staggered g-tensor. 

We have simulated the anisotropic Heisenberg model 
([l]) using the Density Matrix Renormalization Group (DMRG) 
technique 1^]. We have, in general, investigated carefully 
the dependence of the results on the number m of states 
kept in the DMRG calculations. If not stated otherwise, 
we have used open boundary conditions, the finite-system- 
size algorithm and extrapolated the finite-size data to the 
thermodynamic limit, setting J xy = 1 in (|l]). 
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Fig. 1. Comparison between DMRG results (filled circles) for 
the finite-size gap A(L), using m = 50 states, and the Bethe- 
Ansatz results (lines). The gap is given as a function of the 
anisotropy J z in the absence of external magnetic fields. 



1.2 Comparison of DMRG and Bethe-Ansatz results 

In order to verify the accuracy of our approach for the 
excitation energies we use the Bethe-Ansatz equations |1| 
[l9| to evaluate the gap A{L) for finite chains with size L, 

A(L) = E X {L) - Eo(L) , 

for the anisotropic antiferromagnetic Heisenberg chain and 
then compare it with our DMRG simulations. In Fig. [j] 
we plot the finite-size gap obtained by using DMRG ver- 
sus the anisotropy (0 < J z < 1), the lines correspond to 
the exact finite-size gap obtained using Bethe-Ansatz. Al- 
though in this regime the gap in the thermodynamic limit 
is equal to zero, the correspondence between our DMRG 
simulations for the gap for finite system sizes with the 
Bethe-Ansatz results is very accurate. 



2 Homogeneous transversal magnetic field 

Uimin et al proposed, by a mean- field calculation |ic[| , that 
a gap opens for the spin-1/2 anisotropic Heisenberg chain 

L-l 

H = Jxy(Si S* +1 + SfSf +1 ) + J Z S*S* +1 
i=i 

L 

-h x ^Sf (2) 
1=1 

in the presence of a homogeneous longitudinal field h x . 
This mean-field result would then imply that no further 
magnetic anisotropies would be needed to explain the ex- 
perimentally observed spin-gap of Yb4As3. 

In spin-wave theory Uimin et al found Jl0| for the 
model (||) a gap A$wt consisting of two branches: 

Aswt = mm(Ai,A 2 ) , 



Ao 



\ 




(3) 



Up to a certain magnetic field the value of the gap 
is almost linear and above this value the gap starts to 
close itself following a quadratic form, compare Fig. |i| and 
Fig. | 

The physical reasoning for the gap present in the mean- 
field results reviewed above is the following: neglecting the 
SfSf +l coupling term in the Hamiltonian (|2|) it becomes 
identical to the Ising model in a transversal field, which 
has a two-fold degenerate ground-state and a gap. It has 
been argued 116,171, that this reasoning remains valid also 
for (@). 

We compute the energy gap A(L) using the DMRG. 
We calculate the ground state and the lowest excited states 
energies Ei (i ~ 0, 1, 2) and form the following energy 
difference A\ = E% — Eq and A 2 — E 2 — Eq respec- 
tively between the first and the second excited state with 
the ground state energy. Fig. ^ and Fig. |^ show the be- 
havior of Ai(L) with the transverse magnetic field h x for 
J z = 0.75, 0.9706. 

In the thermodynamic limit the first energy differ- 
ence Ai — > in the low magnetic field regime (h x 
2). This behavior is in agreement with the fact that the 
GS is doubly degenerated in this regime when L — > oo. 
Thus to estimate the eventually system-gap in this regime 
we calculate A 2 . A 2 shows a more monotonous behav- 
ior. For fixed h x , A 2 (L) is decreasing with the system-size 
L (A 2 (100) < 0.04). Fig. | and Fig. | show the extra- 
polate gap for the thermodynamic limit with J xy = 1 
and J z — 0.75, 0.9706. We have evaluated also the in- 
duced uniform magnetization M x , which we present in 
Fig. H for J z — 0.25 and J z = 0.9706 (value which corre- 
sponds to the small anisotropies in xy-plane appropriate 
for Yb 4 As 3 §). 

We observe: 

i) The DMRG data presented in Fig. |] and Fig. || (open 
circles) show a phase diagram divided in two regions 
well separated by a critical magnetic field, which de- 
pends on J z and which is around h x —i 2 for J z ~ 1 



in agreement with the isotropic case |20 . 
For h x below the critical magnetic field the system ap- 
pears gapless. For h x above the critical magnetic field, 
a linear gap opens corresponding to the classical fer- 
romagnetic phase, polarized along x-direction. 
As it appears in the Fig. |] and Fig. |^, where we plot 
the spin- wave theory gap (open diamonds), there is a 
substantial disagreement between our DMRG-simula- 
tion and SWT prediction @. 

We find that the induced magnetization M x satu- 
rate for h x larger than the critical field, as illustrated 
in Fig. | for J z = 0.25 and J z = 0.9706 (and the 
staggered magnetization reduces to zero pHl). This ex- 
plains the opening of a gap linear in h x in this phase. 

Dmitriev et al have used scaling arguments fl7]| and 
found for the Hamiltonian (||) that a gap A ~ (h x y 



ii) 



iii) 
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Fig. 2. Finite-size gaps observed for the Heisenberg chain of 
size L in a transversal magnetic field h x with J xy — 1, J z = 
0.75. 




Fig. 3. Finite-size gaps observed for the Heisenberg chain of 
size L in a transversal magnetic field h x with J xy — 1, J z = 
0.9706. 
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Fig. 4. The gap observed for the Heisenberg chain in a 
transversal magnetic field h x with J xy = 1, J z = 0.75. The 
circle symbols come from DMRG simulations (m = 50) ex- 
trapolated to the thermodynamic limit. The diamond symbols 
represent induced gap calculated from spin-wave theory. 
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Fig. 5. The gap observed for the Heisenberg chain in a 
transversal magnetic field h x with J xy = 1, J z = 0.9706. The 
circle symbols come from DMRG simulations (m = 50) extrap- 
olated to the thermodynamic limit. The diamonds represent 
the induced gap calculated from spin-wave theory. 
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Fig. 6. The magnetization M x along x for a anisotropy J z — 
0.25, 0.9706. Note the saturation of M x above a certain critical 
field. 
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opens for small external magnetic field, with an exponent 
v depending on the anisotropy. Our simulations for small 
anisotropy are not supporting these predictions fl7|,|22||. 
But due to the difficulty to resolve accurately Ai using 
DMRG, we cannot exclude completely an eventual small 
gap. Although till now we are not able to give a clear 
answer to the gap-opening question, but at least it seems 
that if there is a gap, the prefactor of the scaling law would 
need to be very small and in any case too small to explain 
experimental results on Yb4As3. 
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3 Staggered magnetic field 
3.1 Model 



Fig. 7. 

fields h: 



Heisenberg case J 



xy 



L-l 



H 



Jr. — J ■ 



L L 

■^Et-w-^E-^- (4) 

i=l i=l 



The above Hamiltonian is not invariant under reflection 
with respect to the mid point of the chain when L is even. 
However in the standard implementation of the DMRG 
algorithm, L is even and the reflection symmetry is used. 
The Hamiltonian (jj) can be easily made invariant under 
reflection by means of a local rotation, given by: 



100 

L 



DMRG results for the the gap for various staggered 



tag 



as a function of chain length L, for h z 



0. 



The staggered Dzyaloshinsky-Moriya interaction, which is 
allowed |j, 7],[nJ in the 4/-compound Yb4As3, leads in 
an external homogeneous magnetic field h z to an effec- 
tive transversal staggered field h x ta9 . In this context the 
anisotropy is not relevant and we can consider the (isotropic) 
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8. Magnetization curves along the z-direction, compare 
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The transformed Hamiltonian reads: 



L-l 



H — J [ Sf Sf +1 (Sf Sf 
=1 



CV QV \ 



L L 

h^Y^S? - h z J2Sf (5) 

i=l i=l 



3.2 Effect of the staggered transversal field 

We start by analyzing the finite-size dependence of the en- 
ergy gap. The Fig. R shows the behavior of the finite size 
system gap versus the staggered magnetic field. One can 
clearly observe (i) that the gap vanishes in the thermody- 
namic limit only for zero h^ ta9 and (ii) that the magnetic 
correlation length is large but finite in the gapped case 
fostag ^ q. data for larger system size L is essentially 
flat for h s x tag > 0. 

In Fig. H we show the magnetization along the z-direc- 
tion as a function of h z , for various h^ tag . Since J xy = 



J x = 1 in ([i]), these results can be compared directly to 
those given in Fig. ^ for fields along x-direction. We note 
that the second-order phase transition to a completely 
magnetized state occurring at h z — 2 in the absence of 
a staggered field is progressively smeared out by h x tag ; 
the transversal field induces quantum fluctuations into the 
magnetized state. These results confirm a similar study 

M 



3.3 Comparison with experiment 

Now using the experimental estimation for exchange cou- 
pling in the isotropic Heisenberg chain J ~ 26 K, we fit 
the experimental data for the gap. The staggered magnetic 
field h x tag is proportional to the experimental magnetic 



field H ext via h s J ag = c sin(0) H e 



c /g± sin(9)h z 



where g± is g-factor for a magnetic-field perpendicular 
to the chain-direction. The g-factors for Yb4As3 are very 
anisotropic, g± has been estimated to be Q g± ss 1.3. 

Our fit yields co/g±sm(8) ~ 0.18. The fitting to ex- 
perimental data depends on cq which has not yet been de- 
termined precisely by experiment. While in the literature 
[ pd| Co — 0.27, we get a good agreement with our DMRG 
simulations and the experimental curve for Co — 0.23 (as- 
suming g ± ~ 1.3)(Fig§. 
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Fig. 9. Fit to experiment with data (m = 70) extrapolated to 
the thermodynamic limit. 



3.4 General case 

We consider now the case where the staggered and the 
uniform magnetic field are not proportional. We want to 
explore the excitation gap of 

L-l L L 

i—1 i i 

_ (6) 

In Fig. Owe present for h z = the gap A as a function 
of the anisotropy 8 and the parameters ao and a\ entering 
in the scaling-law 



A 



«o 



(hi ta9 ) 



(7) 



for the gap. Bosonization predicts jl2],[R| clq w 1.85 and 
a\ = 2/3 for the isotropic case (5 — 0, h z = 0) for small 
staggered magnetic field hf ag . For the isotropic case (8 = 
0) we present in Fig. |ll| the same data as a function of h z . 

We also examined the DMRG-data using (@) includ- 
ing multiplicative logarithmic corrections |1J,[L|: A = 

ao(h s J a9 ) ai \\ogh s J; a9 \ 1 ^ We found essentially the same 
values for the parameters ao and a\ as presented in Fig. [To] 
and Fig. [ll]. We note, that (^) hold only in the asymp- 
totic limit h s ^ ag — * 0, which wc do not examine in the 
present study. The results presented in Fig. [TTJ and Fig. [fl] 
show ao and ai as obtained for overall fits to the gap, for 
y^tag ^ Q g -y^ e b e ii eve this parameter-region to be exper- 
imentally relevant. 



4 Conclusions 

We have studied the anisotropic Heisenberg-chain with 
staggered and uniform transversal and uniform longitudi- 
nal field by DMRG. We found no evidence for a substantial 
gap opening for a homogeneous transversal field, as pre- 
dicted by a mean-field |l0| and a scaling JlTj analysis. We 
found, however, a gap opening for a staggered transversal 
field, in accordance with previous studies [jl2[|l3|| . These 
results lead to the conclusion, that the gap-opening in 
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Fig. 10. Left panel: For h z — the gap versus anisotropy 8 in 
the thermodynamic limit (m = 40) . Right panels: Dependence 
of the parameters ao and a\ as a function of 8. 
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Fig. 11. Left panel: For 5 = the gap versus uniform longi- 
tudinal field h z in the thermodynamic limit (m = 40). Right 
panels: Dependence of the parameters ao and ai as a function 
of 8. 



the the rare-earth, 4/-compound Yb4As3 compound in 
an external magnetic field is attributable to a staggered 
Dzyaloshinsky-Moriya tensor p], p[[7[] ■ 
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